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SUMMARX 



The object of this thesis is to determine the effect of small 
openings on the stress distribution in webs of beams subjected to pure 
bending and to devise a rapid means for predicting the increase in stress 
over that of the condition of no opening. 

Experimental evidence for a family of specimens which covered 
the usual structural practice has been plotted in the RESULTS . A theor- 
etical solution* developed in the APPENDIX , is included and compared with 
the experimental, evidence. From these plots the authors conclude thatt 
the theoretical solution compares quit© favorably rdth actual experiment; 
the stress concentration factor for the edge of trie hole varies directly 
with the ratio of depth of opening to depth cf web and inversely with 
the ratio of length to depth of openingj and that the stress concentration 
factor for the edge of the web varies directly with the ratio of depth 
of opening to depth of web, but the variation with the ratio of length 
to depth of opening is more complex showing a pronounced maximum in the 
vicinity of 1 . 5 * The fact that the ratio of depth of opening to depth 
of web is of little importance until it reaches the value of about 1/3 
is considered excellent verification cf presently accepted practice. 

On the basis of the foregoing, it is considered highly desirable 
that further investigation for specimens having a ratio of length to depth 
of opening in the neighborhood of 1.5 be nude with a view to confirmation. 
It Is further considered desirable that a systematic survey of a similar 
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family be made to determine the effectiveness of reinforcement at open- 
ings and simultaneously cases of complex bendingj i.e., bonding with shear. 
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IHTRODUCTIOH 



In almost all structures openings of son© sort are required, 
either for access of for passing piping, ventilation, or electrical systems 
to the various compartments. It, therefore, behooves the structural 
designer to make proper provision for these requirements without adversely 
affecting the structure. 

Naturally, the circle, because of its simplicity, was the first 

* 

type of opening chosen. Kir sch^ analysed the effect of a circular 

(2) 

opening in a uniform tension field; Tuzi v ' performed the analysis for 

/ 

the circular opening in pure bending; Rowland and Stevenson w/ extended 

both the Kirach and Tuzi analyses for finite width of plate and introduced 

(4) (5) 

the analysis for complex bending. Wolf and Seuber performed analyses 

( 6 ) 

for the elliptical opening in pure bending; Inglis made a similar analysis 

for the elliptical hole in a uniform tension field; Wolf also extended 

his analysis to cover the ellipse in complex bending. Quite recently 

Greenspan^ introduced an analysis for the flat ovaloid in a uniform 

( 3 ) 

tension field. In a forthcoming publication' , Joseph and Brock develop 
a theoretical solution for the flat ovaloid in pure bending. This thesis 
is devoted to the latter problem. The theoretical solution in the APPENDIX 
follows the method proposed in (8). The effect of finite depth is evolved 
by an analysis of experimental data rather than by a rigorous mathematical 
method. 



# Numbers in brackets refer to the Bibliography. 
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It is the purpose of thi3 thesis to correlate data fro® photo- 



elastic experiments with a systematic series of models with varying ratios 
of length to depth of opening and depth of opening to depth of web, sub- 
jected to pure bending, in which the openings are flat ov&loids; i.e., 
openings composed of rectangles cith semicircles mounted at tne long- 
itudinal end3. This correlation includes a check on the mathematical 
solutions developed in (2), (3)> (8), and in the APPENDIX * 



PROCEDURE 



The models were prepared using a nominal thickness of 3/3 
inch ©3 clear Catalin. All models were cut from the same sheet of material, 
and a standard tensile specimen was included. The tensile specimen was 
used only to determine the fringe constant of the particular sheet to 
ensure satisfactory sensitivity as well as to provide means for obtain- 
ing checks through the equations of static equilibrium, if necessary. 

Those models were machined using high-speed steel end mills for all cutting* 

In an effort to minimize the time-edge effect, polymerization, 
a specimen was tested immediately after machining, allowing only suffic- 
ient time for machining stresses tc dissipate. A white light polariscope 
was used for this purpose. The specimen was fitted with end brackets 
and suspended in an A - frame, as shown in Fig. I, whioh subjected the 
entire specimen to pure banding. Suspension was accomplished by vertical 
links fitted at each end with circular pins of drill stock to assure as 
nearly point sup;>ort and loading as was feasible. The A - frame was loca- 
ted in the center of the circular polariscope using a mercury vapor lamp 
with a Wratten #62 Filter. The resulting light was very nearly mono- 
chromatic with a wave length cf 5461 A°. The resulting fringes were 
alternately light and dark giving excellent contrast. 

The specimen was first suspended with brackets and no load* 

The analyzer was rotated to give a dark field. A picture was taken in 
this condition, hereinafter referred to as "no load". The negative 
was then developed to assure proper focus before load was applied. 
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Load was then applied, and a new picture token. The load was previously 
determined which would give a maximum stress of approximately 100 pal., 
well within the elastic range yet sufficiently high to give about seven (7) 
fringes. The analyser remained in the same position as in the nc load 
condition. The negative was then developed to assure satisfactory focus 
and contrast before the specimen was removed for further machining. By 
following tills method, it was possible to complete a run, five (5) specimens 
- all with the same ratio of length to depth of opening, but with varying 
ratios of depth of opening to depth of web, in about four (4) hours which 
minimised the time-edge effect. 
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RESULTS 



The results of the experimental data are incorporated in Figures 
II # III, IV, and V. Figure III is a plot of stress concentration factor 
at edge of hole versus the ratio of length to depth of opening for contours 
of depth of opening to depth of web. Figure II is a plot of stress con- 
centration factor at the edge of the hole versus the ratio of depth of 
hole to depth of web for contours of constant ratios of length to depth 
of opening. On this plot the various experimental data is sliown. Simil- 
arly Figure 17 and 7 show respectively the same plots with stress concen- 
tration factor at the edge of web as the substituted variable. Table I 
is a compilation of experimental data from which the aforementioned curves 
have been plotted. 



TABLE X 



flumjaary of experimental data obtained photoelastic ally 



Ratio of length 
to depth of 
opening 


Ratio of depth 
of opening to 
depth of beam 


Stress Concentration Factor 


ft/b 


b/c 


Bto of hole 


Kdere of beam 




0.1115 


2.067 


1.060 




0.239 


2.027 . 


. 1.033 


1,0 


. 0.364 


2.060 


1.024 




0.485 


2.180 


. 1.040 




0.600 




. 1.100 




0.100 ... . . 


1.430 


. . . 1.045 _ 




... 0.200 


1.440 


1.016 


1.5 


0.300 


1.520 


_ 1.045 




0.400 


1.648 


1.124 




0.500 


1.910 






0.100 


1.320 . 


1.026 




0.200. 


1.322 


1.035 


2.0 





1. 32Q 






0.400 1.413 1.072 

0*522 1 x ^22 1.162 





0.100 


1*160 


1.077 




0.200 


1.176 . . . 


1.026 


2.5 


0.300 


1.300 


1.054 




0.400 




1.087 




kJflCL 


1-235 


UJS 
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DISCUSSION OF RESULTS 



The curveB shown in Figures II and III, and IV and V respectively 
represent the cross-fairing of the data incorporated in Table I. The results 
of Figures II and III are consistent with the well-verified facts that 
a small crack in the longitudinal direction of a beam causes a stress con- 
centration of unity and a small crack in the transverse direction of the 
beam causes an infinite stress concentration. This character iB shown 
again in the slope of the curves of Figures III. It is not surprising, 
either, that the circular opening causes a greater stress concentration 
than any of the other flat ovaloids (no ovaloids transverse to the axis 
of the beam were tested). This result is similar in content to a compar- 
ison of (1) with (2) and of (6) with (4). 

The curves of Figures IV and ? are not so easily explained. 

At the outset there appears to be no rational reason for the apparent 
maximum stress concentration factor at the edge of the beam for a ratio 
of length to depth of opening of 1 . 5 . The almost perfect agreement of 
the experimental results with the theoretical solution for the value of 
the ratio of depth of opening to depth of beam of zero and for all values 
of suoh ratio for the circular opening lend confidence in the remaining 
experimental work. On this basiB the authors feel safe in predicting that 
the maximum augmentation of stress caused by a small opening of this type 
will occur when the ratio of length to depth of opening is in the vicinity 
of 1 . 5 . This same confidence enables the authors to conclude that the 
ratio of depth of opening to depth of beam should lie between 1/4 and l/3 
if no material increase in the stress concentration at the edge of the 
hole is desired. 
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CONCLUSIONS 



1. The stress concentration factor at the edge of the hole varies 
inversely with the ratio of length to depth of opening and directly with 
the ratio of depth of opening to depth of beam. 

2. The stress concentration factor at the edge of the beam varies 
directly with the ratio of depth of opening to depth of beam, rut in some 
complex manner with the ratio of length to depth of hole. 

3. The critical value of line ratio of length to depth of opening 
falls in the vicinity of 1.5. The critical value of the ratio of depth 
of opening to depth of beam lies in the neighborhood of 1/4 to 1/3. The 
critical value is that value at which a small increase in the aforementioned 
ratio is accompanied by a relatively large increase in the stress concen- 
tration factor at the edge of the hole. 

4. The theoretical solution developed in (8) and in the APPENDIX 
is satisfactory for purposes of engineering design. 

5. The photoelastic method is well-suited for inexpensive and rapid 
determination of stress concentration factors. 



RECOMMENDATIONS 



1. Further experimental investigations for values of the ratio 
of length to depth of opening in the vicinity of 1.5 seers highly desir- 
able to fix this critical value. 

2. Further experimental and theoretical investigations in this 
nature to determine the effect of reinforcement of such openings and the 
nature of complex bending) that is* bending with shear, seem desirable. 

3. In the absence of a clearly defined "no load* picture in a photo- 
elastic investigation, it is recommended that small but appreciable stress 
be introduced to clarify the basic pattern. 
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APPENDIX A 



DETAILS OP PROCEDURE 
Description of Models 

Four (4) basic models were prepared. Each was fitted with an 
interior hole l/4 inches in the transverse direction, but the longitud- 
inal extent were respectively, l/4 inches, 3/8 inches, l/2 inches, and 
5/8 inches. The distance between suspension points was maintained at 11 1/8 
inches; the distance from suspension points to the loading points was 
maintained at 5 inches. Each of these models was successively diminished 
in depth in order to obtain the increasing values of the ratio of depth 
of opening to depth of beam. Only in the case of the circular opening 
was it necessary to recut the interior hole. The particulars of the various 
models and the loads to which they were subjected are summarized in Table II. 
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TABLE II 



Summary of Model Particulars 



Hole Dimensions Depth of Beam b/c Load In lbs. a/b 



Ci.UU 

1/4* x 1/4* 


VU %/ X. X-f V 

2.240 


0.1115 


-AM. 

35.35 


1.0 


1/2 1/2 


2.094 


0.239 


35.35 


1.0 


1/2 x 1/2 


1.372 


0.564 


46.00 


1.0 


1/2 x 1/2 


1.030 


0.485 


22,00 


1.0 


1/2 x 1/2 


0.833 . 


. 0.600 


4.00 


1.-0 


1/4 x 3/8 


2.500 


0.100 


85.35 


1.5 


1/4 x 3/3 


1.2500 


0.200 


20.00 


1.5 


1/4 X 3/8 


0.833 


0.300 


9.00 


1.5 


1/4 uc 3/3 


0.625 


0.400 


4.0 


1.5 


1/4 x 3/8 . 


0.500 


0.500 


_ 1.50 ... 


1.5 . 


1/4 x 1/2 


2.500 


0.100 


35.35 


2.0 


1/4 x 1/2 


1.250 


0.200 


20.00 


2,0 


1/4 x 1/2 


0.833 


0.300 


9.00 


2.0 


1/4 x 1/2 


0.625 


0.400 


4.00 


2.0 


_ 1/4 x 1/2 


0.500 


0.500 


1.50 


2*Q__ 


1/4 x 5/8 


2.500 


0.100 


56.90 


2.5 


1/4 x 5/8 


1.250 


0.200 


18.50 


2.5 


1/4 5/8 


0.833 


0.300 


9.50 


2.5 


1/4 x 5/8 


0.625 


0.400 


6.00 


2.5 


1/4 X 5/8 


0.5C0 


0.500 


' 5.00 





* 

\ 



-lft- 



Methods of Analysis 



Standard photoelastic methods were used, supplemented where 
required with direct analytic methods. A transverse section was chosen 
a minimum of two lengths removed from the longitudinal end of the hole. 

* 

The condition of stress at this section should, according to St, Venant 
be removed from the influence of the hole. This section was designated, 
accordingly, "pure bending", see Figs. VI, VII, VIII, and H, and was the 
reference condition for determining the stress concentration factors. 

Examination of the photographs. Figs. X, XI, XII, and XIII, 
revealed that the maximum stress at the edge of the hole occured at its 
vertical axis of symmetry, designated "diameter"} whereas the maximum 
stress at the edge of the beam occured vertically above the longitudinal 
end of the opening, designated "end of hole." It was therefore necessary 
to determine the fringe order existing at the free boundaries in these 
two locations. Since a fringe is the locus of points of constant differ- 
ence between principal stresses, one of which oust be zero at a free boundary, 
no other experimental information was necessary for the evaluation of actual 
stress existing at these free boundaries. Fortunately, the "exact," math- 
ematical, stress analysis for pure bending also reveals that one principal 
stress is also zero everywhere, which fact obviated any additional inform- 
ation for the pure bending condition as well. Thus the desired stress 
concentration factor was merely the ratio of the fringe order actually 
measured at a given point to the fringe order determined at the corres- 
ponding depth in the pure bending condition. Both fringe orders referred 

* Ref. ( 9 ) » p. 31 
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to above must be the actual value determined by the photograph corrected 
for any residual stress in the "no load* condition. 

This correction is shown in Figures VI, VII, VIII, and H as 
the difference between the dotted and solid lines. The case of pore 
bending was obtained by plotting integral fringe orders, solid lines in 
the photographs, at the proper fraction of depth from the bottom. The 
upper and lower halves of the beam were handled separately, since, in pur© 
bending, one half would be in tension and the other in compression, with 
zero stress at the middle surface. The best straight line, guided here 
by the theoretical linear distribution of bending stresses, was passed 
through the experimental points in each half of the beam. It will be 
observed in Figures VI, VII, VIII, and II that in no case did these lines 
pass through zero at the Biddle surface. 8hifts of both halves of these 
lines affected to make the lines pass through this zero point and to inter- 
cept the same value at the edges of the beam. It will be noted from the 
figures that in moat cases this was almost a parallel shift of lines with 
approximately the same values for both halves of the beam. This amount 
of shift, in all but the shallowest models, was of the same character as 
the load stress and less than fringe order 1. This represents the residual 
stress in the specimen, which fro® the density of blackness in the *no 
load* picture, must be very nearly constant over the section. It may 
be noted, however, that in the shallowest beams there is appreciable rosi- 
dual stress, more than a single fringe order and of the same character 
to the load stress. This condition is caused by the bending effect of 
th8 end brackets and admits a more straightforward correction procedure. 
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The Hiaaeter* and "end of hole* curves have this same correc- 
tion factor applied to them in all cases except for the circular opening* 
Figure VI and Z. This procedure was followed since the density of 
blackness on the *no load* picture appeared to be the same throughout 
these specimens except for insignificant surface blemishes. The circular 
openings retained machining stresses for a considerable period of time. 

The fringes representing these machining stresses were concentric rings 
which permitted a plot of fringe order versus fraction of depth from the 
bottom, 3ee Figure VI. These values were then added algebraically to the 
experimentally determined values, previously plotted. The character 
of the machining stress was easily determined by noting the “bulging* 
or *depressing* effect on the fringes in the loaded condition. For example 
when fringes observed in the compression half of the beam were “depressed* 
toward the hole and bent away from the diameter toward the end of the hole, 
thus increasing the fringe order at the diameter of the hole, the machin- 
ing stress was obviously compression and should be subtracted from the 
experimentally determined value to obtain the actual fringe order caused 
by the bending alone. Conversely, if, in the compression half of the 
beam, fringes wore "bulged* away from the edge of the hole, causing dips 
on either side of the diameter, the machining stress was tension and 
should be added to the experimentally determined value to obtain the actual 
fringe order caused by the bending alone. 

A polar plot to obtain the experimental fringe order at the 
diameter of the circular hole wae used in an effort to ascertain this 
peak value with greater accuracy. This is shown in Figure VI. The 
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photograph wna tacked to a drawing board and the center of the hole deter- 
mined. Radial lines were then passed from the center of the circle through 
the intercepts of the fringes at the circumference. The angles were then 
transferred to the plot and the correct value of fringe order was plotted 
radially on this line. Finally, a smooth curve was passed through the 
points and the value of the fringe order dete 'gained on the diameter ex- 
tended. This gave the desired end point of the curve. 
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APPENDIX B 



Sample Calculations 

The data are summarised in entirety in Table I of RESULTS . 
Figures VI, VII, VIII, and IX, referred to in APPEMDIX A. . are shown to 
illustrate the methods of analysis and typical calculations of stress con- 
centration factors. One sample from each series of models is included} 
that is, one for each ratio of length to depth of opening and each ratio 
of depth of opening to depth of beam. 
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APPENDIX C 



Original Data 

The original data is presented in the form of photographs of 
the *no load" and load conditions. Each page of photographs includes all 
the pictures taken for a given series of aodelat Figure X for the ratio 
of length to depth of opening of 1, Figure XI for the ratio of 1.5* Figure 
XII for 2, and Figure XIII for 2.5» 




«A - l 



FIG.X 
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b/ff - 0.5 No Load 



b/C ■ 0.5 Load 





b/C - 0.4 No Load 



b/C - 0.4 Load 





b/C - 0.3 No Load 



b/C » 0.2 No Load 



b/C - 0.1 No Load 



b/C ■ 0.3 Load 




b/C «0.2 Load 
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FIG.ZZIX 




b/C - 0.5 No Load 




b/C » 0.4 No Load 




b/C » 0.5 No Load 




b/C = 0.2 No Load 




b/C =0.5 Load 




b/C - 0.4 Load 




b/C = u.5 Load 




b/C - 0.2 Load 








b/C - 0.5 No Load 




b/C - 0.4 No Load 



a/b = 2.5 



FIG. 2IK 




b/C “ 0.5 Load 




b/C ■ 0.4 Load 




b/C “0.3 Load 




b/C « 0.2 Load 




b/C “0.1 Load 
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APPENDIX D 



K omenc 1 a t u re 



(x,y) = 
< a ,P ) = 
z = 

V = 



5 = 
h = 
X = 

a* 

3 , t, R = 

r*i — 

Jl — 

§.§ = 
9#^ = 
F = 

0 x* °y* a a' °p ~ 

T xy # T ap 

A n» B n» a n» ~ 
M = 

y 0 = 

I = 



SCF ~ 



Cartesian co-ordinates 
Orthogonal curvilinear coordinates 
x / iy 
a / ip 
e lf 

stretch ratio 
e 1 > 

angle in z-plane between tangent to curve 
p = constant and x-axis 
real parameters 
t / jp - >G 

potential functions of complex variable z 
potential function of complex variable 
Airy stress function 

normal stresses: o is normal to surface for vrtiich 

x = oonstant, etc. 

shear stresses 

complex constants 

apolied bending moment 

half-height of opening at x =0 

8 - t / H 

moment of inertia of beam 

bar above, indicates conjugate complex 

stress concentration factor 
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I evelopm.^nt of Theoretical Folihiora 

Consider a function: z= f { £ ) 

which transforms continuously points in the £ ( a , p } pl^ne 
into points in the z(x, y) plane. Here z = x / iy and 

^ = e a+1 P = e v . 

In particular* se? Figs. XV (a) and (b), points on the 
circle of radius e a , with a held constant, defined by a 
position angle § , will transform into points (x,y) on a 

closed cuj’ve defined by the functional form (1). The radial 
line defined by p = constant transforms into a curve orthogonal 
to those for ''hich a = constant. The case of a = 0 is im- 
portant; it defines the unit circle, y , in the K -plane 

and transforms, by the proper selection of the form of (1), 
into the shape of the opening to be studied. A special symbol 
X Kill be assigned to those £- values vhich lie on this 

circle : 

X = e lf} (2) 

so that the conjugate of X : X = e~* ^ = 1/ X. 

The transformation defined by Hcuation (1) describes the 
relation between an Increment dz in the z- plane and an increment 
di: in the ^-plane. The ratio, expret s"d in polar form: 

-- = he ^ (3) 

dv 

3 - This development follows that of (8) and v^s largely guided 
by ir. J. S. Frock. 
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defines h the stretch r 
tangent t > th^ curve 

P 

Article 2.3-2* h and X 



.tic, and X , the angle between 'Che 
_ constant and th« x-axis; see {.l 6 \ 
are obtained by ci if f epentiating ( 1 ): 



» 



tic 

du 





( 4 ) 



Enuatlm {3) and. {■■'<): 

h elX _ „ i~ _ „ r'U ) (5) 

4 * 4 

where the prime indicates differentiation with respect to the 
variable in the parenthesis. 

Thus, the conjugate of (5) in: 

h e"* 1 ^ - z’(" } (6) 

~ a A 

provided (1) has real coefficients and parameters. 
Multiplication of (5) and (6) gives: 



bX z * ( 4 ) z • ( 4 ) 4 4 =z ‘ ( 4 ) 2 ‘ ( 4 ) e c 

The division of (5) "by (6) fives: 

e 2 lX _ zM.4J e 2ip. 

z ‘ ( 4 ) 



(7) 



18) 



In general, in pli-ne problems of elasticity, the tech- 
nique which follows can be applied to any real rational map- 
ping function z = f{ 4 ); that is, where z is the ratio of two 
polynomials in 4 bavin real coefficients. 
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Consider the map-:- inn function due to tree ns pun l7J : 



/, a s / / - d -w 

*■» Y ^ 

S' 



( 9 ) 



i/here s, t, aid 1 „re re^l constants that are determined 
from the proportions of the opening, and z — x / iy and 

: e a+ip 






Then: x 



/ iy a se a • / te~ a • e“^P / Re"*^ a • e - '^? 



= (se a / te“ a )cor p/ Re”^ a cos3p / i(se G - te” a )sin p 
- iRe”^ G sin3 p. 



Equating the real and imaginary parts separately: 
x = (se a / te“ a ) cosp / Re” ;a cos3 p 



( 10 ) 



y = (se a - te” a ) sinp - Re” a sin3 



P 



When a a 0, Equations (10) define the general oval old opening: 



x = (s / t) cos q / R cos 3 



y = (b - t) einp - R PirO p 



(10a) 



Several special cases are of Interest: 



(a) 


t'hen 


zs 

ii 

o 


, an 


ellipse is 


defined. 




(b) 


Then 


Eat 


S3 0 , 


a circle 


of radius n 


is defined. 


(c) 


When 


ct 

II 

O 


and 


F. < 0, a square rlth 


fillets at the 




corners Is 


approximated . 






(d) 


When 


o 

1) 

•p 


and 


R > C, the 


Bnuare is 


rotated 45°. 



In the general case, by the proper choice of s, t, and R, the 
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equations can be made to approximate the cor: on structural 
openinge , e pec Ificailv the one cooposetl of a recta ru le ’1th 
semicircles at each en . 

f- 0 13 

Fron (9): z 1 ( 4 ) “ e - ~ - ^irr (11) 

^ C* 

how h and X can be determined from (?) and (8). '.'hen 
a = 0, corresponding t> the opening, the subs c riot zero 
is used ^nd from (7): 

h2=z'( X),-( X) - (s _ t, ^ _ 2^, 

( 12 ) 

= s 2 + t 2 +9R 2 - (2fct - 6at)ooc2 P - 6RSCOP4? 

In order to obtain strprr.es in plane problems of 
elasticity by using the complex variable method, two potential 
functions must be determined. These functions can be expressed 

uG : 

§<*> = X A n * n / £ a n /, n 

( 1 ?) 

^(z) = Z F n ?. n t Z b n /z n 

where, A n , a n , ? n , and b n are constants , possibly couple::. 

The solution is thus reduced to the evaluation of there con- 
stants. The stress condition at one boundary, infinity, will 
determine A and B ; the stress condition at the other boundary, 
the opening, will determine a and b . Once the two potential 
functions are completely d etermined . the stresses in the plane 
nay be determined by different!, tion. Thus, these functions 
are counterparts of the Airy strife function. In fact, the 
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Airy sire's function correr. ending to those potential functions 
is given by: 

F(x,y) = Real Part of (j>(z) / J^(z)cz] (14) 

The stresses can be found in terms of the two functions: 



a x / a y =2^$'(z) / J'(z)] 



ay - a x / 21 %y =2 [2 $*(s) / J •(*)) 

4 

The usual transformation equations c. ive: 



o + a = 0 +0 

a p x y 



(a p - a a + 2lT ap )e'^ a y - o x + 2iT xy 



(15) 



(16) 



By combining Foliations (3), (15), and (16): 



°Q + °P 



= 0 



[§•(?) / i'(z)j 



u?) 



- 0 + 21 1 =2 r E ( 4 ) $" ( z ) + J 1 ( s )] ---5i e-1 P 

a ap L a 1 ( 4) 



where the functions of z ar c to be ex cress eel in terms of 



£ (a , p ) : 



<f>(z ) = $ (z( £)] = cp ( £) 

$(*) = $£( £)) = <P ( t) 



Differentiation is with respect t j the variable in ,,urentn^- 
s fir : 



4 - Ref. [id} , p. 135 



(19) 



(z ) 

$' U) 



z'«S) 

= hi5i 

7 '(?) 

= .g : i si .„ 



L z ’k>l 



If, at the opening, the normal stress a and the sne^r 

d 

stress t are zero, which is the c se at a fr^e boundary, 
a P 

it can be shown that: 



<p ( x) + -S-i-.-l- V ’(- 1 -) + 5 (-*-) = 0 (20) 

zM^i x x 

The functions <p ( £) and ^ ( £) may be evaluated 
vith the use of the Cauchy Theorems regarding pol°s and. reel- 

S 

dues . 

From [12] , these are : 



Theorem I. If f( is continuous in the closed region 

and anal;/ tic in the interior Tr lth the possible exce >- 
tion of the point £ * 0, where f ( £) VlaS the form 

f( 5) ~ r ^ / / g(£ ) 

and g( £ ) is analytic; then: 



-J- ' /If !.*> d x - . £i - i: . ... 



where y is the unit circle y = 



for | > 1 

ft 
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4 - n«f. \u] . 






Theorem II. Ii f(£ ) In continuous in the cloned region 
\4\> 1 and analytic in the region e.cterior o A , with 
the possible exception of the point £ = co , where f( s ) 
has the form 

f U ) = A 0 ^ A l^ / A 2 s / • • • / A n i 1 / f “ 



Then: 



----- f IL}2<1\ = -f ( 4 ) + A 5+ aJ+ • • . +A^ " 

2 % 1 J r A -4 o i 4. n 

for |4|> 



Consider Equations (13): 



( z ) = 
" ( z ) = 



^ . n-1 

Z. nA z ♦ 

i *1 



* 



Z n( n-1 )A M z n ~'“-* Z 



n 



n(n+l )a n 
Tn +2 



$' (z) * 



T ,n-l ? nb n 

V n V + f - :h*i 



Now let z approach infinity In the x-d irectlon: 

^ ,(z) x = 0 o=\Z "An^jx — 

$''U) M = ^ Z n(n-l)A n z n j 
t'(z) = [" £ nB n z n_1 ] 

X co L 1 n x = 00 



x. — -o 



( 21 ) 



The boundary conditions at x 



= 00 are: 



o 



x ' I y; c i 



- 0 



y 



a- = 0 
xy 



Inserting E uations (21) into the first of Equations (15) 
along vita these boundary conditions gives: 

’ T y= p ^n (x E 1 y^ n ’ 1 / inA n (x - iy) n "5 x ^=o 

This requires that A n = 0 for all values of n except 2. For 
n = 2 

- fJ - 2 C ?iA 2” - 2iA 2 y] 

A = 1 JL 
2 1 31 

Inserting Equations (21) into the second of Equations (15) 
gives: 

__ >c n— 1~1 

^y= 2 \_(x - iy)*2*l«-A 2 / Z nB n (x / iy) J 



V S5 oO 



This requires that B n = 0 far all values of n except 2. For 
n = 2 




Therefore, for pure bending, regardless of the shape of the 
opening, the potential functions are: 



l (z) - qg 2 * % a r/*n 

$(z) = . i| y 2 , £ V z „ 



(22) 



On the boundary of this opening: 



— A 5 — 



<p 



<p ( A ) =^(f + jr 

31 X X' '(Sx+i+iJjn 

X X J 



I 



na 



n 



(r ) = [- "C 1 +t\ + r V ) - Z - - t | 

x L 41 >• • q + t x + ». x ) n+1 -l 



1 z ' ( 



4 , ( a ) =- =-<s X+ 5 + 3,r + £ 

til X X-' 



1 ( B X ♦ 2 + i' ) n 
_ x x 3 



’ ,( X> If* ♦* X + KX -)% Z --- . ? n 

Ui ' <| +tx + R x ) 



-5. , 2 



O '-X3 



>n 



The boundary condition (20) becomes: 

4 . „ 2 



cp ( X ) 



h X + tX^ + Rp! IM^e ^ ^ 
L4 t I + A + 



O 0,0 

p. x )*Z — 

1 (S +tx 
x 



x"!tx 3 ?"^) 



+ ili(| + t X + RX 3 ) 2 + t 5 b 

81 ' (? +tx + R X 5 ) 



Expanding anc.i collecting terms: 



(23) 



\|f (X) - --(e 4 +t 2 +R 2 -st) - -- ( 2 s t + 2Rt - t 2 - 2 He ) 

4 1 81 

- “ ! ( 2 rb - 2 Rt) x 4 + “p 2 X 6 - i: : (S£i-i£«.:-S!) 

Si SI 81 



X 



18 2Ils 



Hu 



£±<E _ Y _ 

31 ' (dtx +Hl) nfi 

A 



n (E_x^_+ t_x£+_R 



+ ^ (2 + t X + Px 3 ) n 



= 0 



At £ =0, <p ( £ ) has the form: 

/ x n . ( \ i’/2st+2Pt— s 2 \,i 2Rs 15a i 

* = *' «> + aj< -5 > + aj • --5 + g‘ 
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r—4 I ,< 



From Theorem I: 



__1__ r 9(V)d* 
2 Til J X-£ 



r 



81 2 
Oi. £ 



iM 



-Jl. T *3>I 

211 J X 



xjax 

x~4 



a - _ i]j{2£t+_2^t_-sf ) 
31 x ^ 



Opa P 



,-ua 



a- 



81 4 * s 2 £ 

for | £ | > 1 



2 Re 



IK 

81 * 4 71 



for ) 4l > 1 



At ^ = oo t (p ( ^ ) has the form: 



_ w „2 v 2 



D „ 



<p( 5) = -- s“ Kf + Z - 



81 



? 



From Theorem II: 



1 

57i 


r = 

^ x-$ 


- 9(4 ) + -- G 2 
81 


F 


for \ ^ l > X 




r 








_ 1__ 
2 rci 


X jp ( X ) <■“ X _ 

J \ _ V 

\- s 


- 49 ( 4 ) + -- 

81 


s 2 « 3 


for i^|>l 




r 








From the 


first of the 


expressions derived from e 


fach theorem: 


9(4) 


= i>i tv f i 


s t_ +2Rt_-e 2 2 Re 


n ♦ -h 


(24) 


31 *- 


F 5 r 


iJ 





From the second of the express lone derived from each theorem; 



iM r 2 J5 2e t+ 2Rt -s 2 2Rs~l 






?<p ‘ « > = |j L>‘ « J + — ■ i 



V 



(25) 



Comparing ( 2 A ) and {25): 



- 4 ? - 



0 



Putting T = t / 



Rt 


1- 


c 


iks r 


'sjg 


^-i 

Ml 





K ■ 



( 26 ) 



But the conjugate of Equation (2>'’) la: 
1 ! 1 ' 

<p + -wiV' (X ) + +< x ) = o 



Z ’ (^ ) 

Using Equation (26), this becomes: 



iMs r b 
ill L ? ; 



•w + T * 2 + R 
4 1 L 2 >6 



* x^l + i!Js 

J 41 



, -v 2T +^x + - T1 - x ) 

'■ x - o • a’—rrrm 

X~ X 4 



+ \|f ( X ) — 0 

4R 



*<X) = --- [--- +TX 2 + F. X 4 - -V- + t X 
L 2>‘ 2 sx 4 -tx -3R 

At ^ = 0, £) has the form: 

* '? > = ^f(;2 2 ) + B («) 



4 * X 



From Theorem I : 



i r <W X)ax 



IMg, R , 

\ * J 



hi 2 £ 



2*t “ 

r 

At Z> = 00 , ^ ( £ ) has the form: 



for jh> 



* 



«>- --- j" 2 + R £ - !5,1::I 2 : J < 0 «,<? ♦ » 

hi L S£ -t£ - "H J 



From Theorem II: 



i _ r*u)*\ 

2 Tti ^ A-£ 



, . . IF G 

'I' ( S) + r-- T 

41 



^ + R 



s£ 4 - 2T - -5 

4 x z 

£ 

s£ - t£ - 3R 



(s+ t£ 2 + R£ 4 >J for \^|> 
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Thei'efore : 

¥ <S> = 5 |-["- £ ? + -5 ~ + 4 - -4-— --s--3£-(s + t^ + -K 



C -zK - 



3 



(27) 



Using the first of Equations (17) and no tin/ that at a free 
boundary a = 0 , o Q = o 

i-c*n »‘in 

p a =0 ^ n=n 



a -0 



zU) z’(^) 



But 



,-x IHsTsA 2 T R~] , - fix 

^ x)=i i it l™* :s + :^J und v ( x )= ’ 

A, A 

<p‘ (X)« Ab'fi Qa - 23| - L^J and 9 • (1) = 



X' 

, ,, x t ?n 

z U) " E • x5 • X" 



and z* ( ^ ) 



1 ;‘P r S . 2 -, d ] 

■ ,-r L'x^ + iX +nX J 

- -^Tx-’-itp.x 5 ] 

in *- * J 

2 4 

= s - t X - 3H x 



Thus ; 

( 



w 

0 )=z [v,s - x -fS-g iv. f-- Tx3 -^ 3 l 

0p i=o 2 L«r • — r ■ « • r.Tj: w> J 



x 



(o ) = "hi j"(s 2 + oHt-2tT-l'R 2 )sin|-(2Rt-st+s 2 )sin3^+Rssin5pl 
P a=0 i t(s 2 +t^ + 9F 2 ) + ( 6Rt - 2st)cos2P - 6Rscos4p J 



(23) 



( G ) is a maxiraura vhen p = 90°. 
P a=0 



Therefore: 



SCF = 



(o R ) 

f_a= 0 _ 

~VElQ~ 



73 "> 



_ 6 

’ yo 



2s^-3Rs-12.F~+2t(4R-t)(l+g) 
(b +t - 3H) 2 
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If vp nov substitute Equations (1), (9), (26), (27), and 
(11) into Equation (l’O, integrate and combine terns, we 
get th® Airy stress function: 



F(a,p) 



= ^ (e 2 +Jt-2Rt+3‘u 2 )e a +(3Rl 



n,2 

A.— 

2 



-2T* -AR-)e 



.[< 



-f3a, -3a 

i te 


-6a , -7a' 

BTe + P. e j 


[ nin p 


t + Rt -8?v 


3g /St 2s? 

+ \-r - --- - 


Apt > - 


3 


6 3 


3 


a -bal 


r -3a 


Pe ” a 


+ Rte 


sin3 p 4-IHse 


— e 

7 



Rs 

o 



-6a 



J ein5 

(29) 
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Evaluation of Parameters 



It Is necessary to determine the values of H, s, and t 
that will ensure that the mapping function ()) will suit the 
actual opening when gC« C. Several conditions may be im- 
posed. on (10a) in order to determine h , s, t such xv : 

(a) s / t / R « r- horizontal dimension = -}a 

(b) s - t / H = h vertical dimension = >b 

(c) Area of opening = (a - b)b / .’,b^R — Jy&x 

(d) Radius of curvature at end = b 

Eauations assoolatinp the desired variables with the 
geometric properties of the opening must no*' be developed to 
provide a systematic procedure for evaluation. To simplify 
this oroblen we introduce the Greenspan notation: 

■V ( P / 0 ) = s 

(30) 

- ( P - 0 ) r= t 

Or: 

s / t = P 

( 30 a) 

s - t s r* 



Equations 


(10a) become 


X s F COSp 


/ * 


1 cos 3p 


y ss f *t si np 


- | 


’< sin 3p 



Imposing the condition of eoual ar°as, T *e obtain: 

A J y dx 53 (a - ’o)b / ib^ for a > b ( 31 ) 

Substituting (10b) into (31), integrating, and simplifying 
gives : 
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n (P. - 3 H 2 ) * (a - L)b / £b 2 * 



(32) 



Reference to any text on differential calculus shove that 
the radius of curvature at any point is given by the ex- 



ore er ion: 



Radius = 



(1 / y ‘ ) J 



,2x3/2 



Now: y* = f ^ and y“ = 

J dx d p dx J 

But , = q cos r - 3R cos 3 B 

dp v v 



ax 



= gy: 



d x 



(33) 



and 

^ t; P sin p / 3R sin 3 p 



Hence: y» = - &-S2L&.Z-25-S2 °-2J 

? sin P / 3B sin 3 P 



( 3 / 4 - ) 



and 



„ _ -Pd_/ 2?R -3R(P / Ocob 4 P / 6K(P-C,)cos 2 P 

y - + 



(P sin p / 3R sin 3 p ) 2 



(35) 



Substituting (^4) and (35) into (33) and simplifying gives: 

RadluB _ ( P 2 Af )/( - 1 .P 2 /3Pn-:V4R/.Hi 2 )cos2 P-3H( 1 >+.- )coe 4 ftpR 2 
-?Ci/ 27 R 2 - 3 R(P/:,.)coBi> ft^6R(?-Ci)co32 p 



?Iov p = 0 for the longitudinal end of the opening. 



The ref ore : 
p^i us 

P / 9R 



The remaining two equations are simply: 
P / n = la 
Ci t R = jVb 



(37) 



(38) 

(39) 
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G one lilo ring 1 only Filiations (32 ), (37 )* ( ~ 2 ) , and (39) as 
a 8 ye ten of equations for the determination of P, 0,, one- R, 
vp n ov choose b = ? and involve the ratio of length to depth 
of opening, u/b . The rye ten of equations become: 

T / R =a/b 
/ R =1 

o 4 ('K» 

Pq - 3 IT = -(a/b - 1) / 1 
? 71 

1L-Z.2LU: = i 

P / 9R 

Any system of more equations than there are unknowns ia in- 
determinate . Therefore, the four possible combinations of 
(4o) taken three at a time were solved for P, q,, and B. for 
the particular case of a/b = 2. This particular case vac 
chosen since it had been used by Greenspan £ ?] as an il- 
lustration and a ready comparison van available. The outline 
of the opening ae defined by there four solutions are plotted 
in Fig. XVI a no compare a with the actual opening; for which 
the fit was desired. 

Because the ratio of depth of opening to depth of bean 
is an important criterion in the determination of the stress 
concentration factor, the depth of opening seemed the most 
logical first choice to be eld constant. Next in importance 
as governing variable was the length of opening. To complete 
the system of equations for solution for P, Q, and R either 
the radius of curvature at the longitudinal end or the area 
must be chosen. A comparison of the fit showed that either 
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viuld be satisfactory, and a comparison of the tvo values of 



stress concentration factor examined vltk Greenspan's value 
as the yardstick shoved that either v mid b ft sat isfactory. 

A good c-iSft could be made for holding the end radius of 



curvature ecu-O to the luii' height of opening since this 
system gave values of stress concentration factor that were 
slightly high - pessimistic. Vhen plotted on a curve similar 
to Fig. Ill, however, the closeness of fit of values calcu- 
lated by holding the area constant vus c lightly superior to 
those values calculated by holding the end radius of curvature 
equal to the half height of opening. The e...s e of calculating 
the values of ?, h, and P * r hen holding the area constant v.-us 
a further consideration. 



Consequently a system of equations consisting of the 
first three of Equations (ko) 1f as chosen: 

P = a/b - H 

Q, = 1 — 1 . ( 4 0 ) 

PQ-OR 2 = i(a/b - 1) / 1 



From uhich: 

R = - ,(a/b/l) / l: /(a/b/1) 2 - 3(4/it- l)(a/b-l) (40a) 



These values of P, Q,, and R are plotted versus a/b in Fig. 

XVII . 

This should not be construed to mean that th»ee values are 
precise; they are merely the best choice that the authors could 
make vith existing information. Somp other geometric property, 
the peripheral length, may warrant consideration. 



- 55 - 



such as 



APPENDIX E 



g. \ b I lQ^rap h y 



(1) Kirsch, G.* "Die Theorie De r Blast! si tatuud Die Bediirfnisao Der 

Featigkaitslehre", Zeitschrift Des Vereinea Deutscher 
Ingenieure, Vol. 42, 1898. 

(2) Tuzi, Z.i "Effect of a Circular Hole on the Streaa Distribution 

of a Beam Under Uniform Bending Moment", Scientific 
Papers of the Tokyo Institute of Physical and Chemical 
Research, Vol. 9, Aug. 20, 1928, pp. 65 - 89. 



(3) Howland, R. C. J., and Stevenson, A. C.x "Biharraonic Analysis in 

a Perforated Strip", Philosophic Transactions of the 
Royal 8ociety of London, Series A, vol. 232, 1933# 

pp. 155 - 222. 

(4) Holf, K.x "Beitrage zur ebenen Elastisitatstheorle", Zeitschrift 

fur technische Phyoifc, vol. 2, no. 8, 1921, p. 209, 
and vol. 3# no. 5 , 1922# pp. 160 - 166. 

(5) Neuber, H . t "Theory of Notch Stresses", David Taylor Model Basin, 

Translation Hr. 74# by A. Raven and J. S. Brock, 
Nov,, 1945 (Distributed by J. I. Edwards, Ann Arbor, 
Mich.) 

(6) Inglis, C. E.x "Stresses in a Plate Due to the Presence of Cracks 

and Sharp Corners", Transactions of the Institution 
of Naval Architects, vol. 55# London 1913# pp* 219 - 
242. 



(7) Greenspan, M.t "Effect of a Small Hole on the Stresses in a Uniformly 

Loaded Plate", Quarterly Applied Mathematics, vol. 2, 
1944# PP. 60 - 71. 

(8) Joseph, J. A. and Brock, J. S.t "The Stresses Around a Small Opening 

in a Beam Subjected to Pure Bending", Preprint of 
Paper No. 50 - APM-3 for presentation at the Conference 
of the Applied Mechanics Division, Purdue University, 
Lafayette, Ind. , June 22 - 24, 1950# of the American 
Society of Mechanical Engineers. 

(9) Timoshenko, S.x "Theory of Elasticity*, McGraw-Hill Book Co., Inc., 

New York, 1934. 



Bibliography (Contra) 



(10) Coker, E. G. and Filon, L. N. G.i *A Treatise on Photoelasticity*, 

Cambridge University Press, London, 1931. 

(11) Sokolnikoff, I. 8.* "Mathematical Theory of Elasticity*, (mimeo- 

graphed lecture notes) , Brown University, Providence, 
Rhode Island, 1941. 



59 - 



DA2B DUB 




































-■ 














































































1 

_ . . 




























\ 





























AUG 3 I B I UDERY 



Thesis 

K144 


12&66 

TTarl 


AUTHOR 

The effect of small holes on the 
T £tpess dis tribut ion-in webs 
subjected to pure bending 


DATE DUE 


BORROWER'S NAME 












J 1 — 



Thesis 12866 

E144 Karl 

The effect of small holes 
on the stress distribution 
in webs subjected to pure 
bending;. 



OR iTSIItl 
r iioui uernaaw 

AW» 

v SMITH 
iil- ; umMW'H 



i.'o 

U. S. Naval Postgraduate School 
Monterey, California 



